Nasopharyngeal angiotensin converting enzyme 2 (ACE2) expression as a risk-factor for SARS-CoV-2 transmission in concurrent hospital associated outbreaks.

BACKGROUND: Widespread human-to-human transmission of the severe acute respiratory syndrome coronavirus two (SARS-CoV-2) stems from a strong affinity for the cellular receptor angiotensin converting enzyme two (ACE2). We investigate the relationship between a patient's nasopharyngeal ACE2 transcription and secondary transmission within a series of concurrent hospital associated SARS-CoV-2 outbreaks in British Columbia, Canada. METHODS: Epidemiological case data from the outbreak investigations was merged with public health laboratory records and viral lineage calls, from whole genome sequencing, to reconstruct the concurrent outbreaks using infection tracing transmission network analysis. ACE2 transcription and RNA viral load were measured by quantitative real-time polymerase chain reaction. The transmission network was resolved to calculate the number of potential secondary cases. Bivariate and multivariable analyses using Poisson and Negative Binomial regression models was performed to estimate the association between ACE2 transcription the number of SARS-CoV-2 secondary cases. RESULTS: The infection tracing transmission network provided n = 76 potential transmission events across n = 103 cases. Bivariate comparisons found that on average ACE2 transcription did not differ between patients and healthcare workers (P = 0.86). High ACE2 transcription was observed in 98.6% of transmission events, either the primary or secondary case had above average ACE2. Multivariable analysis found that the association between ACE2 transcription (log2 fold-change) and the number of secondary transmission events differs between patients and healthcare workers. In health care workers Negative Binomial regression estimated that a one-unit change in ACE2 transcription decreases the number of secondary cases (ß = -0.132 (95%CI: -0.255 to -0.0181) adjusting for RNA viral load. Conversely, in patients a one-unit change in ACE2 transcription increases the number of secondary cases (ß = 0.187 (95% CI: 0.0101 to 0.370) adjusting for RNA viral load. Sensitivity analysis found no significant relationship between ACE2 and secondary transmission in health care workers and confirmed the positive association among patients. CONCLUSION: Our study suggests that ACE2 transcription has a positive association with SARS-CoV-2 secondary transmission in admitted inpatients, but not health care workers in concurrent hospital associated outbreaks, and it should be further investigated as a risk-factor for viral transmission.

Targeted genomic sequencing of avian influenza viruses in wetland sediment from wild bird habitats.

Diverse influenza A viruses (IAVs) circulate in wild birds, including highly pathogenic strains that infect poultry and humans. Consequently, surveillance of IAVs in wild birds is a cornerstone of agricultural biosecurity and pandemic preparedness. Surveillance is traditionally done by testing wild birds directly, but obtaining these specimens is labor intensive, detection rates can be low, and sampling is often biased toward certain avian species. As a result, local incursions of dangerous IAVs are rarely detected before outbreaks begin. Testing environmental specimens from wild bird habitats has been proposed as an alternative surveillance strategy. These specimens are thought to contain diverse IAVs deposited by a broad range of avian hosts, including species that are not typically sampled by surveillance programs. To enable this surveillance strategy, we developed a targeted genomic sequencing method for characterizing IAVs in these challenging environmental specimens. It combines custom hybridization probes, unique molecular index-based library construction, and purpose-built bioinformatic tools, allowing IAV genomic material to be enriched and analyzed with single-fragment resolution. We demonstrated our method on 90 sediment specimens from wetlands around Vancouver, Canada. We recovered 2,312 IAV genome fragments originating from all eight IAV genome segments. Eleven hemagglutinin subtypes and nine neuraminidase subtypes were detected, including H5, the current global surveillance priority. Our results demonstrate that targeted genomic sequencing of environmental specimens from wild bird habitats could become a valuable complement to avian influenza surveillance programs.IMPORTANCEIn this study, we developed genome sequencing tools for characterizing avian influenza viruses in sediment from wild bird habitats. These tools enable an environment-based approach to avian influenza surveillance. This could improve early detection of dangerous strains in local wild birds, allowing poultry producers to better protect their flocks and prevent human exposures to potential pandemic threats. Furthermore, we purposefully developed these methods to contend with viral genomic material that is diluted, fragmented, incomplete, and derived from multiple strains and hosts. These challenges are common to many environmental specimens, making these methods broadly applicable for genomic pathogen surveillance in diverse contexts.

Targeted genomic sequencing with probe capture for discovery and surveillance of coronaviruses in bats.

Public health emergencies like SARS, MERS, and COVID-19 have prioritized surveillance of zoonotic coronaviruses, resulting in extensive genomic characterization of coronavirus diversity in bats. Sequencing viral genomes directly from animal specimens remains a laboratory challenge, however, and most bat coronaviruses have been characterized solely by PCR amplification of small regions from the best-conserved gene. This has resulted in limited phylogenetic resolution and left viral genetic factors relevant to threat assessment undescribed. In this study, we evaluated whether a technique called hybridization probe capture can achieve more extensive genome recovery from surveillance specimens. Using a custom panel of 20,000 probes, we captured and sequenced coronavirus genomic material in 21 swab specimens collected from bats in the Democratic Republic of the Congo. For 15 of these specimens, probe capture recovered more genome sequence than had been previously generated with standard amplicon sequencing protocols, providing a median 6.1-fold improvement (ranging up to 69.1-fold). Probe capture data also identified five novel alpha- and betacoronaviruses in these specimens, and their full genomes were recovered with additional deep sequencing. Based on these experiences, we discuss how probe capture could be effectively operationalized alongside other sequencing technologies for high-throughput, genomics-based discovery and surveillance of bat coronaviruses.

ProbeTools: designing hybridization probes for targeted genomic sequencing of diverse and hypervariable viral taxa.

BACKGROUND: Sequencing viruses in many specimens is hindered by excessive background material from hosts, microbiota, and environmental organisms. Consequently, enrichment of target genomic material is necessary for practical high-throughput viral genome sequencing. Hybridization probes are widely used for enrichment in many fields, but their application to viral sequencing faces a major obstacle: it is difficult to design panels of probe oligo sequences that broadly target many viral taxa due to their rapid evolution, extensive diversity, and genetic hypervariability. To address this challenge, we created ProbeTools, a package of bioinformatic tools for generating effective viral capture panels, and for assessing coverage of target sequences by probe panel designs in silico. In this study, we validated ProbeTools by designing a panel of 3600 probes for subtyping the hypervariable haemagglutinin (HA) and neuraminidase (NA) genome segments of avian-origin influenza A viruses (AIVs). Using in silico assessment of AIV reference sequences and in vitro capture on egg-cultured viral isolates, we demonstrated effective performance by our custom AIV panel and ProbeTools' suitability for challenging viral probe design applications. RESULTS: Based on ProbeTool's in silico analysis, our panel provided broadly inclusive coverage of 14,772 HA and 11,967 NA reference sequences. For each reference sequence, we calculated the percentage of nucleotide positions covered by our panel in silico; 90% of HA and NA references sequences had at least 90.8 and 95.1% of their nucleotide positions covered respectively. We also observed effective in vitro capture on a representative collection of 23 egg-cultured AIVs that included isolates from wild birds, poultry, and humans and representatives from all HA and NA subtypes. Forty-two of forty-six HA and NA segments had over 98.3% of their nucleotide positions significantly enriched by our custom panel. These in vitro results were further used to validate ProbeTools' in silico coverage assessment algorithm; 89.2% of in silico predictions were concordant with in vitro results. CONCLUSIONS: ProbeTools generated an effective panel for subtyping AIVs that can be deployed for genomic surveillance, outbreak prevention, and pandemic preparedness. Effective probe design against hypervariable AIV targets also validated ProbeTools' design and coverage assessment algorithms, demonstrating their suitability for other challenging viral capture applications.

Mutations in emerging variant of concern lineages disrupt genomic sequencing of SARS-CoV-2 clinical specimens.

Mutations in emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages can interfere with laboratory methods used to generate viral genome sequences for public health surveillance. We identified 20 mutations that are widespread in variant of concern lineages and affect widely used sequencing protocols by the ARTIC network and Freed et al. Three of these mutations disrupted sequencing of P.1 lineage specimens during a recent outbreak in British Columbia, Canada. We provide laboratory validation of protocol modifications that restored sequencing performance. The study findings indicate that genomic sequencing protocols require immediate updating to address emerging mutations. This work also suggests that routine monitoring and protocol updates will be necessary as SARS-CoV-2 continues to evolve. The bioinformatic and laboratory approaches used here provide guidance for this kind of assay maintenance.

The contrasting role of nasopharyngeal angiotensin converting enzyme 2 (ACE2) transcription in SARS-CoV-2 infection: A cross-sectional study of people tested for COVID-19 in British Columbia, Canada.

BACKGROUND: Angiotensin converting enzyme 2 (ACE2) protein serves as the host receptor for SARS-CoV-2, with a critical role in viral infection. We aim to understand population level variation of nasopharyngeal ACE2 transcription in people tested for COVID-19 and the relationship between ACE2 transcription and SARS-CoV-2 viral load, while adjusting for expression of: (i) the complementary protease, Transmembrane serine protease 2 (TMPRSS2), (ii) soluble ACE2, (iii) age, and (iv) biological sex. The ACE2 gene was targeted to measure expression of transmembrane and soluble transcripts. METHODS: A cross-sectional study of n = 424 "participants" aged 1-104 years referred for COVID-19 testing was performed in British Columbia, Canada. Patients who tested positive for COVID-19 were matched by age and biological sex to patients who tested negative. Viral load and host gene expression were assessed by quantitative reverse-transcriptase polymerase chain reaction. Bivariate analysis and multiple linear regression were performed to understand the role of nasopharyngeal ACE2 expression in SARS-CoV-2 infection. FINDINGS: Analysis showed no association between age and nasopharyngeal ACE2 transcription in those who tested negative for COVID-19 (P = 0•092). Mean relative transcription of transmembrane (P = 0•00012) and soluble (P<0•0001) ACE2 isoforms, as well as TMPRSS2 (P<0•0001) was higher in COVID-19-negative participants than COVID--19 positive ones, yielding a negative correlation between targeted host gene expression and positive COVID-19 diagnosis. In bivariate analysis of COVID-19-positive participants, transcription of transmembrane ACE2 positively correlated with SARS-CoV-2 viral RNA load (B = 0•49, R2=0•14, P<0•0001), transcription of soluble ACE2 negatively correlated (B= -0•85, R2= 0•26, P<0•0001), and no correlation was found with TMPRSS2 transcription (B= -0•042, R2=<0•10, P = 0•69). Multivariable analysis showed that the greatest viral RNA loads were observed in participants with high transmembrane ACE2 transcription (Β= 0•89, 95%CI: [0•59 to 1•18]), while transcription of the soluble isoform appears to protect against high viral RNA load in the upper respiratory tract (Β= -0•099, 95%CI: [-0•18 to -0•022]). INTERPRETATION: Nasopharyngeal ACE2 transcription plays a dual, contrasting role in SARS-CoV-2 infection of the upper respiratory tract. Transcription of the transmembrane ACE2 isoform positively correlates, while transcription of the soluble isoform negatively correlates with viral RNA load after adjusting for age, biological sex, and transcription of TMPRSS2. FUNDING: This project (COV-55) was funded by Genome British Columbia as part of their COVID-19 rapid response initiative.

Assessing oligonucleotide designs from early lab developed PCR diagnostic tests for SARS-CoV-2 using the PCR_strainer pipeline.

INTRODUCTION: During the first month of the SARS-CoV-2 outbreak, rapid development of PCR-based diagnostic tests became a global priority so that timely diagnosis, isolation, and contact tracing could minimize the advancing pandemic surge. Designing these tests for broad, long-term detection was complicated by limited information about the novel virus' genome sequence and how it might mutate during global spread and adaptation to humans. METHODS: We assessed eight widely adopted lab developed PCR tests for SARS-CoV-2 against 15,001 SARS-CoV-2 genome sequences. Using a custom bioinformatic pipeline called PCR_strainer, we identified all mismatches and sequence variants in genome locations targeted by 15 sets of primer/probe oligonucleotides from these assays. RESULTS: For 12 out of 15 primer/probe sets, over 98 % of SARS-CoV-2 genomes had no mismatches. Two primer/probe sets contained a single mismatch in the reverse primer that was present in over 99 % of genomes. One primer/probe set targeted a location with extensive polymorphisms with 23 sequence observed variants at the forward primer location. One of these variants, which contains three nucleotide mismatches, arose in February as part of the emergence of a viral clade and was present in 18.8 % of the genomes we analyzed. DISCUSSION: Most early PCR diagnostic tests for SARS-CoV-2 remain inclusive of circulating viral diversity, but three assays with extensive mismatches highlight assay design challenges for novel pathogens and provide valuable lessons for PCR assay design during future outbreaks. Our bioinformatics pipeline is also presented as a useful general-purpose tool for assessing PCR diagnostics assays against circulating strains.

The SOS Response Master Regulator LexA Regulates the Gene Transfer Agent of Rhodobacter capsulatus and Represses Transcription of the Signal Transduction Protein CckA.

UNLABELLED: The gene transfer agent of Rhodobacter capsulatus (RcGTA) is a genetic exchange element that combines central aspects of bacteriophage-mediated transduction and natural transformation. RcGTA particles resemble a small double-stranded DNA bacteriophage, package random ∼4-kb fragments of the producing cell genome, and are released from a subpopulation (<1%) of cells in a stationary-phase culture. RcGTA particles deliver this DNA to surrounding R. capsulatus cells, and the DNA is integrated into the recipient genome though a process that requires homologs of natural transformation genes and RecA-mediated homologous recombination. Here, we report the identification of the LexA repressor, the master regulator of the SOS response in many bacteria, as a regulator of RcGTA activity. Deletion of the lexA gene resulted in the abolition of detectable RcGTA production and an ∼10-fold reduction in recipient capability. A search for SOS box sequences in the R. capsulatus genome sequence identified a number of putative binding sites located 5' of typical SOS response coding sequences and also 5' of the RcGTA regulatory gene cckA, which encodes a hybrid histidine kinase homolog. Expression of cckA was increased >5-fold in the lexA mutant, and a lexA cckA double mutant was found to have the same phenotype as a ΔcckA single mutant in terms of RcGTA production. The data indicate that LexA is required for RcGTA production and maximal recipient capability and that the RcGTA-deficient phenotype of the lexA mutant is largely due to the overexpression of cckA. IMPORTANCE: This work describes an unusual phenotype of a lexA mutant of the alphaproteobacterium Rhodobacter capsulatus in respect to the phage transduction-like genetic exchange carried out by the R. capsulatus gene transfer agent (RcGTA). Instead of the expected SOS response characteristic of prophage induction, this lexA mutation not only abolishes the production of RcGTA particles but also impairs the ability of cells to receive RcGTA-borne genes. The data show that, despite an apparent evolutionary relationship to lambdoid phages, the regulation of RcGTA gene expression differs radically.